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ABSTRACT 

We investigate a stationary pair production cascade in the outer magnetosphere of 
a spinning neutron star. The charge depletion due to global flows of charged particles, 
causes a large electric field along the magnetic field lines. Migratory electrons and/or 
positrons are accelerated by this field to radiate curvature gamma-rays, some of which 
collide with the X-rays to materialize as pairs in the gap. The replenished charges 
partially screen the electric field, which is self-consistently solved, together with the 
distribution functions of particles and gamma-rays. By solving the Vlasov equations 
describing this pair production cascade, we demonstrate the existence of a stationary 
gap in the outer magnetosphere of PSR B1055-52 for a wide range of current densities 
flowing in the accelerator: From sub to super Goldreich-Julian values. However, we 
find that the expected GeV spectrum becomes very soft if the current density exceeds 
the Goldreich-Julian value. We also demonstrate that the GeV spectrum softens with 
decreasing magnetic inclination and with increasing distance to this pulsar. We thus 
conclude that a sub-Goldreich- Julian current, a large magnetic inclination, and a small 
distance (500 pc, say) are plausible to account for EGRET observations. Furthermore, 
it is found that the TeV flux due to inverse Compton scatterings of infrared photons 
whose spectrum is inferred from the Rayleigh- Jeans side of the soft blackbody spectrum 
is much less than the observational upper limit. 

Subject headings: gamma-rays: observations - gamma-rays: theory - - magnetic fields 
- pulsars: individual (B1055-52) - X-rays: galaxies 
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1. Introduction 

The EGRET experiment on the Compton Gamma 
Ray Observatory has detected pulsed signals from 
seven rotation-powered pulsars (for Crab, Nolan et 
al. 1993, Fierro et al. 1998; for Vela, Kanbach et 
al. 1994, Fierro et al. 1998; for Geminga, Mayer- 
Hasselwander et al. 1994, Fierro et al. 1998; for PSR 
B1706-44, Thompson et al. 1996; for PSR B1046-58, 
Kaspi at al. 2000; for PSR B1055-52, Thompson et 
al. 1999; for PSR B1951+32, Ramanamurthy et al. 
1995). The modulation of the 7-ray light curves at 
GeV energies testifies to the production of 7-ray radi- 
ation in the pulsar magnetospheres either at the polar 
cap (Harding, Tademaru, & Esposito 1978; Daugh- 
erty & Harding 1982, 1996; Michel 1991; Sturner, 
Dermer, & Michel 1995; Shibata, Miyazaki, & Taka- 
hara 1998; Harding & Muslimov 1998; Zhang & Hard- 
ing 2000), or at the vacuum gaps in the outer mag- 
netosphere (e.g., Cheng, Ho, & Ruderman 1986a,b, 
hereafter CHRa,b). Both of these pictures have had 
some success in reproducing global properties of the 
observed emission. However, there is an important 
difference between these two models: An polar gap ac- 
celerator releases very little angular momenta, while 
outer gap one could radiate them efficiently. More 
specifically, the total angular momentum loss rate 
must equal the energy loss rate divided by the angu- 
lar velocity of the star, implying an average location 
of energy loss on the light cylinder (Cohen & Treves 
1972; Holloway 1977; Shibata 1995), of which distance 
from the rotation axis is given by 

where il denotes the angular frequency of the neutron 
star, and c the speed of light. 

On these grounds, the purpose here is to explore 
a little further into a model of an outer gap, which 
is located near the light cylinder. If the outer mag- 
netosphere is filled with a plasma so that the space 
charge density is equal to the Goldrcich-Julian den- 
sity, pgj = — QB z /(2ttc), then the field-aligned elec- 
tric field vanishes, where B z is the component of the 
magnetic field along the rotational axis. However, the 
depletion of charge in the Goldreich and Julian model 
in a region where it could not be resupplied, may 
cause a vacuum region to develop. Holloway (1973) 
pointed out the possibility that a region which lacks 
plasma is formed around the surface on which pej 
changes its sign. Subsequently, CHRa,b developed 



a version of an outer magnetospheric 7-ray emission 
zone in which acceleration in the Holloway gaps above 
the null surface, where B z vanishes, brought particles 
to large Lorentz factors (~ 10 7,5 ). 

After CHRa,b, many outer-gap models have been 
proposed (Chiang & Romani 1992, 1994; Romani 
& Yadigaroglu 1995; Romani 1996; Zhang & Cheng 
1997; Cheng, Ruderman & Zhang 2000). Their basic 
idea is as follows: A deviation of the charge density 
from pgj results in an electric field along the magnetic 
field, B. If this electric field becomes strong enough to 
accelerate e ± pairs to ultrarelativistic energies, they 
could radiate 7-rays tangential to the curved B field 
lines there. The curvature 7-rays may be converted 
into e 1 ^ pairs via two photon collisions. In order to 
keep a steady current flow and the charge density pgj 
in the regions outside the gap, the gap will grow un- 
til it is large enough and the electric field is strong 
enough to maintain a sufficient supply of charges to 
the rest of the open field line region. If the gap ends 
in a region where pgj 7^ and E B 7^ 0, charges from 
the surrounding region may flow in through the end 
(or the boundary) . If both boundaries are located on 
the null surface, pairs produced in the gap will re- 
place the charge deficiency inside the gap, and finally 
the gap will be filled up. However, if a vacuum gap 
extends to the light cylinder, the charged particles 
created in the gap should escape from the magneto- 
sphere, so the gap would not be quenched. Hence, 
stable outer gaps (if they exist) are those from the 
null surface to the light cylinder along the last-open 
field lines. In an outer gap, the inner boundary lies 
near the intersection of the null surface where pc.T 
vanishes and the boundary of the closed field lines of 
the star on which the magnetospheric current does 
not exits. 

Recently, solving the Vlasov equations that de- 
scribe a stationary pair production cascade, Hirotani 
and Shibata (2001a,b, Papers VII and VIII) eluci- 
dated important characteristics of the particle ac- 
celeration zones in the outer magnetosphere. They 
demonstrated that a stationary gap does not ex- 
tend between the null surface and the light cylin- 
der. Rather, their width (W) along the field lines 
is adjusted so that a single particle emits copious 7- 
rays one of which materialize as a pair in the gap 
on average. The resultant W becomes about 5% of 
wlc m the case of the Crab pulsar. The produced 
pairs in the gap do not completely cancel the accel- 
eration field, because the particles do not accumulate 
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at the gap boundaries. Outside of the gap, the par- 
ticles flows away from the gap along the magnetic 
field lines as a part of the global current flow pat- 
tern. In another word, the remained, small-amplitude 
electric field along the magnetic field lines outside of 
the gap prevents the gap from quenching. Moreover, 
they demonstrated that the gap position is not fixed 
as considered in previous outer gap models. Their 
position shifts as the magnetospheric current flowing 
the gap changes. That is, if there is no particle injec- 
tion across either of the boundaries, the gap is located 
around the null surface, as demonstrated in their ear- 
lier papers (Hirotani & Shibata 1999a, b,c, Papers I, 
II, III; Hirotani 2000a, b,c, Papers IV,V,VI; see also 
Bcskin et al. 1992 and Hirotani & Okamoto 1998 for 
a pair-production avalanche in a black-hole magne- 
tosphere). However, the gap position shifts towards 
the light cylinder (or the star surface) if the injected 
particle flux across the inner (or the outer) bound- 
ary approaches the Goldreich- Julian value. In other 
words, the accelerator can be located at any place in 
the magnetosphere; their position is primarily con- 
strained by the magnetospheric current. 

In Papers VII and VIII, they examined the outer 
gaps formed in the magnetospheres of young or mil- 
lisecond pulsars, assuming a sub- Goldreich- Julian par- 
ticle injection across the boundaries. In this paper, 
we apply the same method to a middle-aged pulsar 
B1055-52 and demonstrate that a stationary gap is 
maintained even if the injected particle flux exceeds 
the Goldreich-Julian value, and that the possibility of 
such a large particle flux injection is ruled out if we 
compare the predicted 7-ray spectrum with observa- 
tions. 

In the next section, we formulate the basic equa- 
tions describing the stationary pair production cas- 
cade in the outer magnetosphere. In § 3, we ap- 
ply the method to a middle-aged pulsar, B 1055-52 
and reveal that large magnetic inclination and a sub- 
Goldreich-Julian current are plausible for the outer- 
magnetospheric accelerator for this pulsar. In the 
final section, we compare the results with previous 
works. 

2. Basic Equations 

We first consider the continuity equations of parti- 
cles in § 2.1, the 7-ray Boltzmann equations in § [2.2] , 
and the Poisson equation for the electrostatic poten- 
tial in 
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2.1. Particle Continuity Equations 

Under the mono-energetic approximation, we sim- 
ply assume that the electrostatic and the curvature- 
radiation-reaction forces cancel each other in the 
Boltzmann equations of particles. Then the spa- 
tial number density of positrons (e + 's) and electrons 
(e~'s), N + (s) and 2V_(s), obey the following continu- 
ity equations: 



dN ± 



where 



Q(x) 



7 [V p +G++Vp- 



_G_ 



(2) 



(3) 



Here, G+(x, e 7 ) and G-[x, e 7 ) refer to the distribu- 
tion functions of outwardly and inwardly propagating 
7-ray photons, respectively, having energy m e c 2 e 7 ; v+ 
(or V- ) designates the velocity of the center of mass 
of e + 's (or e~'s). The pair production rate for an 
outwardly propagating (or inwardly propagating) 7- 
ray photon to materialize as a pair per unit time is 
expressed by ?7 P + (or %>-). For charge definiteness, 
we consider that a positive electric field arises in the 
gap. In this case, positrons (or electrons) are migrat- 
ing outwardly (or inwardly). 

Since gyration will be canceled out in v± , we obtain 



v± = ±ccos$— + rQsin#e <! >, (4) 



where 



1 d 



r sin 9 < 



(5) 



is the azimuthal unit vector, $ = sin _1 (ri7 sin#/c) is 
the projection angle of the three-dimensional particle 
velocity onto the poloidal plane, and B = \B\. Impos- 
ing the stationary condition d t + fld^ = 0, we obtain 
from equations (@) and (0) 



± V I ccos$^iV ± J = 



de 1 [r]p + G + + ry p _G_ 



(6) 

In this paper, we estimate $ at the gap center such 
that 

_ x r cnt fisin6» cn t 



$ = sin 



(7) 



where r cnt is the distance of the gap center from the 
star center, and 9 cn t is its colatitude angle. Utilizing 
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V • B = 0, and assuming that the toroidal bending is 
negligible in the sense that \B$\ <C B, we obtain 



-G- 



ds \ B J ccos$ J 7 p 

(8) 

The pair production redistribution functions ry p ± 
can be defined as 



f? P ±(e 7 ) = C 1 - Mc> 



de x ^^cr p (e 7 ,e x ,fi c ), (9) 



where cr p is the pair-production cross section and 
cos _1 ^ c refers to the collision angle between the 7- 
rays and the X-rays with energy m e c 2 e x . X-ray num- 
ber density between dimcnsionless energies e x and 
e x + de x , is integrated over e x from the threshold en- 
ergy e th = 2(1 - ^c) - ^^ 1 to infinity. 

To evaluate /x c between the surface X-rays and the 
7-rays, we must consider 7-ray's toroidal momenta 
due to the aberration of light. It should be noted 
that a 7-ray photon propagates in the instantaneous 
direction of the particle motion at the time of the 
emission. Therefore, evaluating the 7-ray toroidal ve- 
locity at the gap center, we obtain 



/'c 



cos $ cos / 



poh 



(10) 



where 6 po i is the collision angle between the surface 
X-rays and the curvature 7-rays, that is, the angle be- 
tween the two vectors (r cn t,0 C nt) and (B r ,B ) at the 
gap center. Therefore, the collision approaches head- 
on (or tail-on) for inwardly (or outwardly) propagat- 
ing 7-rays, as the gap shifts towards the star. 

2.2. Boltzmann Equations for Gamma-rays 

Unlike the charged particles, 7-rays do not propa- 
gate along the magnetic field line at each point, be- 
cause they preserve the directional information where 
they were emitted. However, to avoid complications, 
we simply assume that the outwardly (or inwardly) 
propagating 7-rays dilate (or constrict) at the same 
rate with the magnetic field; this assumption gives a 
good estimate when W <C h7lc holds. Under this 
assumption, we obtain (Paper VI) 



ds\B 



ccos$ 



G± 



Vc 



ccos<I> 



where (e.g., Rybicki, Lightman 1979) 



V3e 2 r 1 



F — 



N±, (11) 



(12) 



1 3 hcT 3 
° m c c 2 4-7T Rq ' 

POO 

F(u) = u I Kn(t)dt; 



(13) 
(14) 



i?c is the curvature radius of the magnetic field lines 
and K 5 / 3 is the modified Bessel function of 5/3 order. 
The effect of the broad spectrum of curvature 7-rays 



is represented by the factor F(e 7 /e c ) in equation ([121). 
r refers to the particle Lorentz factor. 

2.3. Poisson Equation 

The real charge density e(N + — N-) appears in the 
Poisson equation for the non-corotational potential 
Neglecting relativistic effects, we can write down the 
Poisson equation as follows: 



V^* = 4tt 



e(N+- iV_) + 



rtB z 



2ttc 



(15) 



where e designates the magnitude of the charge on an 
electron. 

2.4. One-dimensional Analysis 

As described at the end of § 3 in Paper VII, it is 
convenient to introduce the typical Debey scale length 
c/wp, 



Une 2 nB cnt 
m c 2nce 



(16) 



where i? cn t represents the magnetic field strength at 
the gap center. The dimcnsionless coordinate variable 
then becomes 

By using such dimensionless quantities, and by as- 
suming that the transfield thickness, D±, of the gap 
is greater than W, we can rewrite the Poisson equa- 
tion ( |l5| ) into the following one-dimensional form 
(e.g., § 2 in Michel 1974): 



F # 



(18) 



dE \\ 



A 2 



m + ^ ds) 



B, 



where = e*(s)/(m c c 2 ) and Aj_ = (uj p /c)D±; 

the particle densities are normalized by the Goldreich- 
Julian value as 



MO 



2nce N± 



(20) 
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We evaluate B z /B at each point along the last-open 
field line, by using the Newtonian dipole field. 

Let us introduce the following dimensionless 7-ray 
densities in the dimensionless energy interval between 
(3i-i and 0c 

2-Kce pP* 

9±i0=7^- rfe 7 G±( S ,e 7 ). (21) 



In this paper, we set 0q = 10 2 , which corresponds to 
the lowest 7-ray energy, 51.1 MeV. We divide the 7- 
ray spectra into 11 energy bins and put 0i = 10 2 ' 5 , 

02 = 10 3 , 03 = 10 3 - 25 , 04 = 10 3 - 5 , 05 = 10 3 - 75 , As = 

10 4 , 7 = 10 4 - 25 , 0s = 10 4 - 5 , 9 = 10 4 - 75 , 1O = 

10 5 , 0n = 10 5 25 . We can now rewrite the continuity 
quation (|J) of particles into 
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dn± _ 

d£ B cos $ ^ 

s i=\ 



[fb+'g+W + rb-W-iOl ( 22 ) 



where the magnetic field strength, B, is evaluated at 
each £. The dimensionless redistribution functions 
rj p± are evaluated at the central energy in each bin as 



%± 



1 

— ?7 P ± 

LUp 



0i-i + 0i 



(23) 



A combination of equations (|22|) gives the current 
conservation law, 



jtot = "-+(0 + "--(0 = constant for £. 



(24) 



If 

Jtot = 1 holds, the current density per unit mag- 
netic flux tube equals the Goldreich- Julian value, 
fi/(2?r). 

The Boltzmann equations ([ll[) for the 7-rays are 
integrated over e 7 between dimensionless energies 
0i-i and 0i to become 

igl = ± (In B) 9 i T H:9± ± TT^W (25) 
d£ y± d£ K ,y± cos$ y± S cnt cos$ ' v ; 

where i = 1, 2, • • •, m (to = 9) and 

V3e 2 r rft/^ 



7 /c 



Lu p hR c 



du / Ks(t)dt (26) 



is dimensionless. 

Equating the electric force e\d^f/dx\ and the radi- 
ation reaction force, we obtain the saturated Lorentz 
factor at each point as 



/3i? c 



V 2e 



ds 



1 



1/4 



(27) 



we compute the curvature radius R c at each point 
for a Newtonian dipole magnetic field. In the case of 
middle-aged pulsars such as B1055-52, their less dense 
X-ray fields lead to a large pair-production mean free 
path for a 7-ray photon to materialize as a pair. To 
ensure that the Lorentz factor do not exceed the max- 
imum attainable limit, ^(£2)1 we compute T with 



1 
f 



1 



1 



V> 2 (6 



(28) 



2.5. Boundary Conditions 



We now consider the boundary conditions to solve 
the Vlasov equations @, @, (j22j) , and @. At 
the inner (starward) boundary (£ = £1), we impose 
(Paper VI) 

£||(£i) = 0, (29) 

= 0, (30) 

gU&)=° (< = l,2,-..,9). (31) 

It is noteworthy that condition (|2^) is consistent with 
the stability condition at the plasma- vacuum interface 
if the electrically supported magnetospheric plasma is 
completely-charge-separated, i.e., if the plasma cloud 
at ^ < £1 is composed of electrons alone (Krause- 
Polstorff & Michel 1985a,b; Smith et al. 2001). We 
assume that the Goldreich- Julian plasma gap bound- 
ary is stable with E\\ = on the boundary, £ = £1. 

Since positrons may flow into the gap at £ = £1 as 
a part of the global current pattern in the magneto- 
sphere, we denote the positronic current per unit flux 
tube at £ = £1 as 



«+(60 = h, 

which yields (eq. J24|) 

= .itot -h- 



(32) 



(33) 



At the outer boundary (£ = £2), we impose 

£||(6) = 0, (34) 

g i _fo) = (i = 1,2,- • -,9), (35) 

n-(6)=j 2 . (36) 

The current density created in the gap per unit flux 
tube can be expressed as 



Jgap = Jtot -Oi -J2- 



(37) 
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We adopt j gap , ji, and j% as the free parameters 
We have totally 2m + 6 boundary conditions (|2 



(36) for 2m + 4 unknown functions \E', Eh, n±, g± 
(i = 1,2, ■ ■ -,m), where m — 11. Thus two extra 
boundary conditions must be compensated by making 
the positions of the boundaries £i and £2 be free. The 
two free boundaries appear because En = is imposed 
at both the boundaries and because j gap is externally 
imposed. In other words, the gap boundaries (£1 and 
£2) shift, if j\ and/or j'2 varies. 

It is worth mentioning that the gap width, W, is 
related with the X-ray field density as follows (for 
details, see Papers V, VII, and VIII): 



W = 



A 



P Jgap 

jtot ' 



(38) 



where A p refers to the pair production optical depth, 
and Nj the expectation value of the 7-ray photons 
emitted by a single particle while it runs the gap. 
Equation (|8|) is automatically satisfied by the sta- 
tionary Vlasov equations. The relation is useful when 
we interpret the variation of W. For example, for 
middle-aged pulsars, the X-ray field is much less 
dense compared with the magnctosphcric X-ray field 
of young pulsars. We can therefore expect that the 
gap is more extended (i.e., W is greater) for middle- 
aged pulsars than the young pulsars, because the for- 
mer's larger A p leads to a larger W, compared with 
the latter. We examine such features more accurately 
in the next section, by analyzing the Vlasov equations 
numerically under the boundary conditions (B9h-(36). 



3. Application to PSR B 1055-52 
3.1. X-ray and Infrared Field 

Combining ROSAT and ASCA data, Greiveldinger 
et al. (1996) reported that the X-ray spectrum con- 
sists of two components: a soft blackbody with kT s = 
68 eV and A s = 0.78 A (d/0.5) 2 and a hard blackbody 
with kT h = 320 eV and A h = 2.5 x 10- 5 A*(d/0.5) 2 , 
where d refers to the distance in kpc. Here, A s and A^ 
indicate the observed emission erea of soft and hard 
blackbody components, respectively; A* expresses the 
area of the whole neutron star surface. 

Mignani, Caraveo, and Bignami (1997) claimed 
that the blue/UV flux is consistent with an extrapo- 
lation of the soft, thermal part of the ROSAT X-ray 
spectrum. We thus compute the infrared flux from 
the Rayleigh- Jeans side of the soft X-ray spectrum 



with kT s = 68 eV and A s = 0.78A(d/0.5) 2 . This 
IR photon field is, in fact, so weak that the IC drag 
acting on a particle is negligibly small compared with 
the curvature drag, as the TeV/GeV flux ratio (<C 1) 
indicates in the following subsections. Therefore, the 
Vlasov equations can be solved without the IR field 
(i.e., the pair-production rate can be computed only 
from the X-ray field). Once En(£), n±(£), and <?± l (£) 
are solved, we can passively compute the TeV spec- 
trum. 

3.2. Dependence on Injected Currents 

Let us now substitute the X-ray field into equa- 
tion (|^) and solve the Vlasov equations by the method 
described in § |^. The explicit GeV and TeV spectra 
are computed by the method described in § 4 in Paper 
VII. 

3.2.1. Equal Current Injection 

We first consider the case when j\ = fa . In figure 
we present the acceleration field, 



ds 



E n 



(39) 



for the four cases: j\ = ji = 0, 0.005, 0.05, and 
2.0, which are indicated by the solid, dashed, dash- 
dotted, and dotted curves, respectively. In the third 
case for instance, the positronic and electronic cur- 
rents equally flow into the gap across the inner and 
outer boundaries at the rate 0.05r2/27r per unit mag- 
netic flux tube (i.e., 5% of the typical Goldreich- 
Julian value). For all the four cases, we set a\ = 60°. 
Since the solution forms a 'brim' to disappear (fig. 2 
in Hirotani & Okamoto 1998) if j gap exceeds typically 
several percent, we adopt j gap = 0.01 as the repre- 
sentative value. The abscissa designates the distance 
along the last-open field line and covers the range 
from the neutron star surface (s — 0) to the position 
where the disance equals s — 0.35ti7lc = 3.29 x 10 6 m. 

It follows from figure |l| that the gap is located 
around the null surface (where B z vanishes), which 
is indicated by the downarrow. For ot\ = 60°, 
n = 31.9 rad s~\ the null surface is located at s = 
1.59 x 10 6 m from the star surface along the last- 
open field line. The conclusion that the gap is located 
around the null surface for j% = j'2, is consistent with 
what was obtained analytically in § 2.4 in Paper VII. 
We can also understand from figure that W, and 
hence —d^/ds decreases with increasing ji = ji- This 
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is because not only the produced particles in the gap 
(i.e., jgap) but also the injected ones (i. e., J1+J2) con- 
tribute to the 7-ray emission. Because of the efficient 
7-ray emission due to the injected particles, a sta- 
tionary pair-production cascade is maintained with a 
small pair-production optical depth W/X p , and hence 
a small W. It is noteworthy that W does not lin- 
early decrease with j mp /j tot, by virtue of the 'nega- 
tive feedback effect' due to the N~ x factor; that is, 
the decrease in iV 7 due to the reduced W partially 
cancel the original decrease in W. It may be note- 
worthy that iV 7 oc WT oc W 3 / 2 holds (for details, see 
eqs. [8], [11], and [31] in Paper V, along with fig. 6 in 
Hirotani k Okamoto 1998). 

We present the 7-ray spectrum in figure the 
curves correspond to the same cases as those in fig- 
ure [lj In each case, outwardly propagating 7-ray 
flux is depicted, because it is greater than the in- 
wardly propagating one. The observed fluxes are in- 
dicated by open circles and squares. It follows that 
the spectrum becomes soft as the injected current 
densities, j\ = 32, increase. This can be easily un- 
derstood because —d^/ds decreases with increasing 
ji = ji ( n g- 0)- Even though solutions exist for su- 
per Goldreich- Julian currents (e.g., the dotted line in 
fig. ||), the predicted spectra become too soft to match 
the observations. On these grounds, we can conclude 
that a sub- Goldreich- Julian current density is prefer- 
able for us to account for the observed spectrum. 

3.2.2. Current Injection across Outer Boundary 

We next consider the case when the gap is shifted 
towards the star surface by virture of the particle in- 
jection across the outer boundary. In figure |], we 
present —d^/ds for vanishing ji, the solid, dashed, 
and dotted curves represent the cases of {31,32)= 
(0,0), (0,0.25), and (0,0.391), respectively. For all 
the three cases, we set a\ — 60° and j gap = 0.01. 
The abscissa covers the range from s = to s = 
0.35za7lc- There is no solution above 32 > 0.391, be- 
cause —d^/ds distribution forms a brim at the inner 
boundary at this critical value. 

It follows from the figure that the gap position 
shifts inwards as 32 increases. For example, the gap 
center position shifts from 0.225?zjlc for 32 = (solid 
curve) to 0.049ro LC and 0.0084ro LC for j 2 = 0.25 
(dashed) and 32 — 0.391 (dotted), respectively, from 
the star surface along the last-open field line. In ad- 
dition, W decreases significantly as the gap shifts in- 
wards. This is because both the factors jgap/jtot and 
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Fig. 1. — Distribution of the acceleration field, 
-W/ds, for PSR B1055-52 for a ; = 60°. The 
solid, dashed, dash-dotted, and dotted curves 
correspond to the cases of (i ga p Ji J2)=(0. 01,0,0), 
(0.01,0.005,0.005), (0.01,0.05,0.05), and (0.01,2,2), re- 
spectively. The downarrow indicates the position of 
the null surface, where B, vanishes. 
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Fig. 2. — Expected pulsed 7-ray spectra from 
PSR B1055-52 for a-, = 60° and j gap = 0.01. The 
four curves correspond to the same cases as presented 
in fig. 1. In each case, the larger flux for either the 
inwardly or the outwardly propagating 7-rays is de- 
picted. 
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Ap decreases in equation (|38|). For details, see § 6.1 
in Paper VII. 

Let us now consider the Lorentz factor. As the 
maximum value in the gap, we obtain 3.0 x 10 7 for 
32 = (solid curve), whereas 8.7 x 10 6 for ji = 0.391 
(dotted one). That is, in the latter case, T reduces 
significantly even though the maximum of —d^/ds is 
greater than the former. This is because the parti- 
cles' motion becomes unsaturated due to the reduced 
W. Therefore, the mono-enegetic approximation is 
no longer valid for the dotted curve in figure |^. 

The emitted spectra are presented in figure [|. Be- 
cause the Lorentz factor decreases as ji increases, the 
spectrum becomes soft as the gap shifts inwards. For 
the dotted curve, the spectrum would be further soft- 
ened from figure ||, if we considered the unsaturated 
effect of particle motion. We can therefore conclude 
that the gap should not be located well inside of the 
null surface for ct\ « 60°, so that the predicted GeV 
spectrum may not be significanly inconsistent with 
observations. 

3.2.3. Current Injection across Inner Boundary 

Let us consider the case when the gap shifts to- 
wards the light cylinder, as a result of the particle 
injection across the inner boundary. In figure [5], we 
present —dfy/ds for vanishing j 2 ] the solid, dashed, 
dash-dotted, and dotted curves represent the cases 
of C?iJ 2 )= (0,0), (0.25,0), (0.5,0), and (0.585,0), re- 
spectively. For all the four cases, we set a; = 60° 
and jgap = 0.01. The abscissa covers the range from 
s = to s = 1.3n7Lc- There is no solution above 
ji > 0.585, because — d^>/ds distribution forms a 
brim at the outer boundary at this critical value. 

It follows from the figure that the gap shifts out- 
wards as ji increases. For example, the gap center 
is located at s = s cnt = (si + s 2 )/2 = 0.225h7lc 
for ji — (solid curve), while it shifts to the posi- 
tion s = 0.526g7lc and s = 0.854h7lc for j\ = 0.25 
(dashed) and ji = 0.5 (dash-dotted), respectively. 
In addition, W increases as the gap shifts outwards. 
This is because the diluted X-ray field at the outer 
part of the gap increases A p (see eq. Q). The in- 
crease of W with increasing j%, however, does not 
mean that —d^/ds increases with j\. This is because 



the small |pgj| fx \B Z 



at larger distance from 



the star results in a small \d(—d^/ds)/ds\ in the Pois- 
son equation (|l5|). 

The emitted spectra are presented in figure ^. 
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Fig. 3.— Distribution of -dV/ds for PSR B1055-52 
for ai = 60°. The solid, dashed, and dotted curves 
correspond to the cases of (j gap ,:/i j2) = (0.01,0,0), 
(0.01,0,0.25), and (0.01,0,0.391), respectively. 
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Fig. 4. — Expected pulsed 7-ray spectra from 
PSR B1055-52 for an = 60°, j gap = 0.01, and ji = 0. 
The four curves correspond to the same cases as pre- 
sented in fig. 3. In each case, the larger flux for either 
the inwardly or the outwardly propagating 7-rays is 
depicted. If the curve is thick (or thin), it means that 
the outwardly (or inwardly) propagating 7-ray flux is 
greater than the other. 
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The curvature spectrum becomes the hardest for 
(ji J2)=(0,0) of the four cases, because —d^/ds be- 
comes large by virtue of the strong magnetic field at 
relatively small distance from the star. It follows from 
the solid and dashed curves in figures || and || that 
the spectrum does not soften very rapidly when j\ 
increases (fig. |) compared with the case when fa in- 
creases (fig. |). This is because W increases with 
increasing ji (fig. ||), while it significantly decreases 
with increasing j 2 (fig. ||). 

Let us summarize the main points that have been 
made in §§ 3.2.1-3.2.3. The curvature spectrum be- 
comes the hardest when ji — j 2 = 0. In this case, 
the gap is located close to the null surface, where B z 
vanishes. 

3.3. Dependence on Inclination 

It has been revealed that the 7-ray spectrum be- 
comes hardest when the injected current is small (say, 
when ji ~ ]2 ~ 0) and that the EGRET spectrum 
cannot be explained for cti = 60° (as the solid line 
111 figs. |, f§ or § indicates). Since the 7-ray energies 
are predicted to increase with increasing a; in Pa- 
per V, we examine in this section whether the EGRET 
pulsed flux around 6 GeV can be explained if consider 
a larger a\. In figures ^| and ||, we present the accel- 
eration fields and the expected spectra for «j = 80°, 
75°, 60°, and 45° as the dash-dotted, dashed, solid, 
and dotted lines. The current densities are chosen 
as jgap = 0.01, ji = 32 = 0. In each case, the out- 
wardly propagating 7-ray flux is depicted, because it 
is greater than the inwardly propagating one. 

It follows from figure || that a large inclination an- 
gle (75° or greater) is preferable to explain the 7-ray 
spectrum around 6 GeV. The reasons are fourfold (see 
also § 6.2 in Paper VII): 

• The distance of the null surface from the star de- 
creases with increasing a{. 

• The magnetic field strength, and hence the Goldreich- 
Julian charge density, QB Z / (2nc), in the gap increases 
as the distance from the star decreases. 

• It follows from the Poisson equation (|l5| ) that the 
derivative of —<N>/ds (i.e., ~d 2l & /ds 2 ) increases with 
increasing B z (i.e., with decreasing distance from the 
star). 

• By virtue of this increasing derivative, the accel- 
eration field, —d^/ds, at the gap center increases. 
We may notice here that W does not decrease very 
rapidly with increasing ai, because of the 'negative' 
feedback effect due to the N 1 factor in equation (p8) . 
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Fig. 5.— Distribution of -d^/ds for PSR B1055- 
52 for ai = 60°. The solid, dashed, dash- 
dotted, and dotted curves correspond to the cases of 
(Jgap Jij 2 )=(0.01,0,0), (0.01,0.25,0), (0.01,0.5,0), and 
(0.01,0.585,0), respectively. 
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Fig. 6. — Expected pulsed 7-ray spectra from 
PSR B1055-52 for a\ = 60°, j gap = 0.01, and j 2 = 0. 
The solid, dashed, dash-dotted, and dotted curves 
correspond to the cases of j\ — 0, 0.25, 0.5, and 
0.5847, respectively. In each case, the larger flux for 
either the inwardly or the outwardly propagating 7- 
rays is depicted. If the curve is thick (or thin), it 
means that the outwardly (or inwardly) propagating 
7-ray flux is greater than the other. 
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Fig. 7.— Distribution of -d^/ds for PSR B1055-52 
for ai = 80° (dash-dotted), 75° (dashed), 60° (solid), 
and 45° (dotted). The current densities are fixed as 
Jgap = 0.01 and ji = j-2 = 0; therefore, the peak of 
each curve corresponds to the null surface for each 
inclination. 
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As a result of this increased —d^/ds, the spectrum 
becomes hard for a large a\. 

It is worth noting that the Tev flux is always neg- 
ligibly small compared with the observational upper 
limit. This is because the infrared flux, which is de- 
duced from the Rayleigh- Jeans side of the soft surface 
blackbody spectrum, declines sharply as v 2 at small 
frequency, hv <C kT s . 

3.4. Dependence on Distance 

We have assumed that the distance is 0.5 kpc so 
that the observed area of the soft blackbody emission, 
A s = 0.78A*(d/0.5kpc) 2 , may be less than A*. In 
this section, we consider the case when the distance is 
1.5 kpc, as indicated by Taylor, Manchester and Lyne 
(1993), and compare the results with d = 0.5 kpc case. 

In figure 0, we present the resultant spectra for 



at = 80°, j g 



0.01, and ji = j% = 0. The solid (or 



1 9 10 10 10 11 10 12 10 13 



dashed) line represents the 7-ray spectrum when the 
distance is 0.5 kpc (or 1.5 kpc). The normalization 
is adjusted by the fluxes below 0.5 GeV; D± = 2.3W 
(or 7.7W) is assumed for d = 0.5 kpc (or 1.5 kpc). It 
is clear from figure || that the GeV spectrum becomes 
hard if we assume a smaller distance, d. 

The conclusion that a small distance to this pul- 
sar is preferable, is derived, because the less dense 
X-ray field for a smaller d results in a larger A p 
in equation (38), and hence in a larger W and the 
acceleration field, —d^>/ds. However, the solution 
does not vary so much if the X-ray density increases 
(1.5kpc/0.5kpc) 2 = 9 times. This is due to the 'neg- 
ative feedback effect' caused by the TV" 1 factor in 
equation (^8|); that is, the increase of W is substan- 
tially canceled by the increase of iV~, which is propor- 
tional to the product of the curvature-radiation rate 
of a particle per unit length and W. In other words, 
the solution exists for a wide range of parameters and 
does not change very much if the variation of the X- 
ray field density is within one order of magnitude. 



photon enerqy [eV] 

Fig. 8. — Expected pulsed 7-ray spectra from 
PSR B1055-52 for 45° (dotted), a { = 60° (solid), 75° 
(dashed), and 80° (dash-dotted). The current den- 
sities are fixed as j gap = 0.01 and j\ = j 2 = 0. In 
all the four cases, outwardly propagating fluxes are 
greater than (or dominate) the inwardly propagating 
ones. 



4. Discussion 

In summary, we have developed a one-dimensional 
model for an outer-gap accelerator in the magne- 
tosphere of a rotation-powered pulsar. Solving the 
Vlasov equations that describe a stationary pair- 
production cascade, we revealed that the accelerator 
shifts towards the star surface (or the light cylinder) 
as the particle injection rate across the outer (or in- 
ner) boundary approaches the Goldreich-Julian value. 
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Applying this theory to a middle-aged pulsar, B1055- 
52, we find that stationary solutions exist for a wide 
range of parameters such as the injected currents, 
magnetic inclination, and the distance to the pulsar. 
Comparing the expected spectrum with ERGET ob- 
servations, we conclude that a sub Goldreich- Julian 
current density, a large magnetic inclination, and a 
small distance (500 pc, say) are preferable for this 
pulsar. 

4.1. Magnetospheric Current Distribution 

As demonstrated in § 6.3 in Paper VII, the spin- 
down luminosity of a rotation-powered pulsar be- 
comes (see § 6.3 in Paper VII) 



r 2 

Jtot /active 



4, ,2 



(40) 



where fi m is the neutron star's magnetic dipole mo- 
ment; /active represents the ratio of magnetic fluxes 
along which currents are flowing and those that are 
open. ForPSRB1055-52,fi 4 /x^/c 3 = 10 3464 ergs s -1 . 
If all the open field lines are active (i.e., /active = 1), 
jtot ~ 1 is required so that the observed spin-down 
luminosity 10 3448 ergs s _1 may be realized. If a small 
fraction of open field lines are active (i.e., /active <C 1), 
a super-Goldrcich- Julian current (i.e., jtot 1) is re- 
quired. 

However, our outer-gap model requires (and, in 
fact, previous outer-gap models assume) j to t <C 1. 
Thus the deficient current must flow along the open 
field lines that are not threading the accelerator. 

4.2. Comparison with Previous Works 

Let us compare the present methods and results 
with Paper V. In the present paper, Eu, N±(s), and 
G±(s, e 7 ) were solved from the Vlasov equations for a 
non-vacuum gap, while in Paper V only E\\ field was 
solved from the Poisson equation for a vacuum gap, 
with the aid of the 'closure condition' (38), where 
jgap = jtot was assumed. Since the condition ( |38| ) is 
automatically satisfied by the Vlasov equations, the 
electrodynamics solved in this paper (e.g., W, Eu) are 
essentially the same with Paper V, provided that the 
gap is nearly vacuum (i.e., jtot *C 1). For example, 
the peak energy of curvature radiation is predicted to 
be about 1 GeV for a>i — 45° in Paper V, which is 
comparable with the vF v peak of the dotted line in 
figure |8|. It also follows from figure || that the peak 
energy does not strongly depend on this conclu- 



sion is consistent with what was obtained for middle- 
aged pulsars in Paper V (see table 3). In the present 
paper, we can more accurately predict the 7-ray spec- 
trum, which is divided into 11 energy bins, compared 
with the gray approximation adopted in Paper V. For 
example, the 7-ray spectrum is predicted to decline 
very sharply above 10 GeV (fig. ||) even for a large a\ 
(~ 75°, say). 

We finally compare the present work with the nu- 
merical simulation performed by Smith et al. (2001), 
who examined equilibrium charge distributions in the 
magnetosphere of an aligned rotator. In the first 
place, they demonstrated that the Goldreich- Julian 
charge distribution is unstable and collapses to form 
a polar dome containing plasma of one charge and an 
equatorial belt containing plasma of the other sign: 
E • B = inside both of them. These are separated 
by a vaccum gap in which E • B 7^ holds. To apply 
the theory developed in the present paper to the case 
when the Goldreich-Julian charge distribution breaks 
down, we must modify the Poisson equation (|l^) (or 
eqs. |l^] and [eq:basic-2]) so that the electric field 
caused by the charges in the dome and the belt may 
be taken into account. 

In the second place, Smith et al. (2001) demon- 
strated that the creation of electron-positron pairs in 
the gap between the dome and the belt reduces the 
value of E ■ B in the gap so that it turns off at last. 
This is because the charges are attracted to regions of 
the same sign, following the closed field lines, and in- 
creasingly filling the magnetosphere. Instead, as con- 
sidered in the present paper, if pairs are created in 
the open field line region, we may expect global flows 
of charged particles and hence a stationary pair pro- 
duction cascade in the magnetosphere of a spinning 
neutron star. If we extend the work by Smith et al 
(2001) into oblique rotators and apply the present 
method to the gap, we may construct a stationary 
magnetospheric model with pair production cascade 
under the existence of global currents. There is room 
for further investigation on this issue. 

One of the authors (K. H.) wishes to express his 
gratitude to Drs. Y. Saito and A. K. Harding for valu- 
able advice. He also thanks the Astronomical Data 
Analysis Center of National Astronomical Observa- 
tory, Japan for the use of workstations. 
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Fig. 9. — Expected pulsed 7-ray spectra from 
PSR B1055-52 for a; = 80°, j gap = 0.01, and j 1 = 
]2 = 0. The solid (or dashed) line represents the 7-ray 
spectrum when the distance is 0.5 kpc (or 1.5 kpc). 
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